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Abstract. We investigate 2-colour QCD with 2 flavours of Wilson fermion at nonzero temperature T and 
quark chemical potential fj,. From temperature scans at flxed fj, we find that the critical temperature for 
the superfiuid to normal transition depends only very weakly on above the onset chemical potential, 
while the deconfinement crossover temperature is clearly decreasing with fj,. We find indications of a region 
of superfiuid but deconfined matter at high /i and intermediate T. The static quark potential determined 
from the Wilson loop is found to exhibit a 'string tension' that increases at large n in the 'deconfined' 
region. The electric (longitudinal) gluon propagator in Landau gauge becomes strongly screened with 
increasing temperature and chemical potential. The magnetic (transverse) gluon shows little sensitivity to 
temperature, and exhibits a mild enhancement at intermediate fi before becoming suppressed at large fi. 
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PACS. 11.15Ha Lattice gauge theory 
12.38Mh Quark-gluon plasma 



12.38Aw Lattice QCD calculations - 21.65Qr Quark matter 



1 Introduction 



The properties and phase diagram of quantum chrome- 
dynamics (QCD) at large baryon density remain largely 
unknown, despite substantial theoretical efforts. One im- 
portant reason for this is the failure of traditional Monte 
Carlo methods for lattice QCD at nonzero density, due to 
the infamous sign problem. While progress has been made 
in the region of high temperature T and moderate baryon 
chemical potential /is, the region of low T and large fiB 
remains inaccessible to Monte Carlo simulations. Q 

There are, however, QCD-like theories which are not 
affected by the sign problem, at least for an even num- 
ber of flavours Nf, among them QCD with gauge groups 
SU(2) (QC2D) [5i,6,7 , or G2 [8], with nonzero isospin den- 
sity [S] or with adjoint fermions [TU]. These theories allow 
first-principles lattice simulations and may hence be used 
as benchmarks for other methods that are not encum- 
bered by the sign problem, but which may involve uncon- 
trolled approximations. Such methods may include model 
studies such as NJL and quark-meson models (possibly 
augmented by a Polyakov loop potential); effective theo- 
ries valid for example for heavy quarks or at high density; 
or functional methods such as the functional renormali- 
sation group (FRG), Dyson-Schwinger equations (DSEs) 



^ Some progress has recently been made using complex 
Langevin [T][2] and other methods [3l|4], but neither have as 
yet been shown to work for QCD. 



or n-particle irreducible (nPI) methods, which rely on as- 
sumptions about higher-order vertices. 

Among these QCD-like theories, QC2D is the simplest, 
both mathematically and computationally, and hence pro- 
vides the greatest scope for quantitative tests of other ap- 
proaches. It has already been used to test analytic contin- 
uation methods at high temperature and low density [11] 
[T^ and to compare model studies with lattice results (see 
P^fHlfTS] for some recent examples). 

In a recent paper [7] the first steps were taken towards 
mapping out the phase diagram of QC2D in the (/i, T) 
plane, by identifying at least 3 distinct regions charac- 
terised by vanishing or nonvanishing diquark condensate 
(qq) and Polyakov loop (L): 

— a vacuum or hadronic phase at low T and /i, charac- 
terised by vanishing quark number density nq and with 
{qq) = 0, (i) « 0; 

— a superfiuid, confined (quarkyonic) region at low T and 
intermediate to large fi, with (L) w 0, {qq) ^ 0; 

— a deconfined quark-gluon plasma at high T, with {L) ^ 
Oi (qq) = 0; and possibly 

— a deconfined, superfiuid region at large fi and interme- 
diate T, with (L) ^ 0, {qq) ^ 0. 

Rough estimates were presented for the phase boundaries 
between these regions. In the present paper, we study the 
phase boundaries in more detail and attempt to pinpoint 
their location. 

We will also attempt to cast further light on the confin- 
ing properties of the theory at low temperature and the 
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nature of the putative deconfinenient transition at high 
density, by computing the static quark potential in the 
low-temperature region. 

In contrast to quantities which may not be directly 
comparable between theories, the effects of the medium 
on low order Green functions in QC2D may provide a re- 
liable guideline to full QCD. Quark and gluon correlation 
functions are of great interest, as the theory can be fully 
expressed in terms of these. Propagators play a predomi- 
nant role, in particular in continuum descriptions, and in 
some cases their behaviour suffices to shed light on the 
critical physics of the phase diagram, e.g. the deconfine- 
ment transition [T6l[T7llT8lfT9] . In this paper, we will study 
how the gluon propagator responds to both temperature 
and quark chemical potential. 

In Sec. [D we set out the details of our lattice simula- 
tions, including the action, parameters and lattice volumes 
used. Then, in Sec. [3] we study the superfluid to normal 
and deconfinement transition by performing a tempera- 
ture scan at 3 different values of the chemical potential. 
The response of the static quark potential to fi is investi- 
gated in Sec.Sl while in Sec.[S]results for the gluon propa- 
gator are reported. Preliminary results for the gluon prop- 
agator have been reported in [20l[2T| . and for the static 
quark potential in |22) . 



2 Simulation details 

We use a standard Wilson gauge action with two flavours 
of unimproved Wilson fermion, with the addition of a di- 
quark source term to lift the low-lying eigenmodes and 
allow a controlled study of diquark condensation effects. 
Further details about the action and the simulation method 
can be found in 5,6,7]. The results obtained will depend 
on the diquark source j; in the end the j ^ limit must 
be taken to obtain 'physical' resultsH 

We use the same parameters as in [Bl[7], namely (3 — 
1.9, K = 0.168, corresponding to a lattice spacing a = 
0.178(6)fm and a pion mass arriTr = 0.645(8). The lightest 
baryon, the scalar diquark, is degenerate with the pion in 
the vacuum, and at zero temperature we therefore expect 
an onset transition to a superfluid phase at m = m,r/2. 
This has been corroborated in previous simulations [SJIHl 

In addition to the ensembles used and described in , 
we have generated gauge configurations on 16'^ x N^- lat- 
tices with Nr =4-20, in order to study in detail the ther- 
mal transitions at an = 0.35,0.4,0.5 and 0.6. The details 
of these ensembles are given in table [TJ For most tem- 
peratures, two diquark sources ja = 0.02 and 0.04 have 
been used, enabling us to perform a linear extrapolation 
to the j = limit. In the region of the superfluid to normal 
transition, where a linear extrapolation is known to be in- 



In cases where model studies could be carried out with 
j 7^ one might also compare directly results for nonzero j; 
however, most other studies will not contain any explicit di- 
quark source term, so the j -t- limit is crucial. 
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Fig. 1. Diquark condensate (qq) as a function of temperature 
T for chemical potential — 0.35, 0.4, 0.5, 0.6 (top to bot- 
tom). The circles are data extrapolated to j = using a linear 
Ansatz for ja < 0.04; the shaded circles denote the results of 
a linear extrapolation using j = 0.02, 0.03 only. 



valid, two additional j-values have been added to allow 
for a controlled extrapolation. 



3 Phase transitions 



3.1 Superfluid to normal transition 



Figure [T] shows the order parameter for superfluidity, the 
diquark condensate 



{qq) = (V''*^C75T2V' - rcj5T2rn , 



(1) 



as a function of the temperature T, for fia — 0.35, 0.4, 0.5 
and 0.6. Also shown are the results of a linear extrapola- 
tion to j = 0. We can clearly observe a transition from a 
superfluid phase, characterised by (qq) ^ 0, at low tem- 
perature, to a normal phase with (qq) = at high temper- 
ature, with a transition in the region 0.08 < Ta < 0.12 for 
all three values of /x. In order to pinpoint the transition, 
we have performed simulations with 4 different j-values 
in the transition region, and used 3 different functional 
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Table 1. Number of trajectories for different temperatures, T — l/{aNr), chemical potentials fj, and diquark sources j. All 
trajectories have average length 0.5, and the spatial lattice size is Ns = 16 in all cases. 



forms for tlie j — >■ extrapolation, 

Imear: 
power law: 
constant + power: 



{qq) =A + Bj, (2) 
(qq) - Br , (3) 
(qq) =A + Br ■ (4) 



The results of these extrapolations are summarised in ta- 
ble [5] We find that the linear form is clearly disfavoured; 
however, the pure power law does not give a good fit ei- 
ther in most cases, while the constant -I- power Ansatz is 
unstable and tends to give a negative intercept at j — 0. 
Clearly, more work is needed to obtain good control over 
the J — extrapolation. 

With our current data at only a single volume we are 
not in a position to determine the order of the transition, 
although it would be expected to be a second order phase 
transition in the universality class of the 3-dimensional XY 
model. This could be tested by attempting a universal fit 
of the data at j 7^ in the transition region to a scaling 
form given by the appropriate critical exponents: note that 
at fixed ^, the diquark source j would be the magnetic field 
variable in this scaling function. This would also provide 
an alternative method for determining Tg. 

Because of the uncertainties in the j ^ extrapola- 
tion in the critical region, we have estimated the critical 
temperatures Ts for the superfluid to normal transition by 
determining the inflection points for (qq) at ja = 0.02 and 
0.04, and extrapolated the resulting values to j = using 
a linear Ansatz. The results are given in table [31 We see 
that Tg is remarkably constant over the whole range of /x- 
values considered. The indications are that the transition 
happens at a somewhat lower temperature at /xa = 0.35, 
but this point is already very close to the onset from vac- 
uum to superfluid at T = 0, /XqC = = 0.32, sug- 
gesting that Ts{iJ,) rises very rapidly from zero ai ^ — Ho 
before suddenly flattening off. 



Fit 






A 


B 


a 




© 


0.4 


9 


0.0036(1) 


0.695(3) 




26 




0.4 


10 


0.0056(1) 


0.711(4) 




102 




0.4 


11 


0.0096(1) 


0.646(4) 




9.8 




0.5 


9 


0.0037(2) 


0.829(4) 




51 




0.5 


10 


0.0095(2) 


0.790(6) 




54 




0.5 


11 


0.0155(2) 


0.712(6) 




45 




0.6 


9 


0.0034(2) 


0.966(4) 




46 




0.6 


10 


0.0095(2) 


0.976(6) 




48 




0.6 


11 


0.0203(2) 


0.849(6) 




180 




0.4 


9 




0.490(8) 


0.852(5) 


7.4 




0.4 


10 




0.438(7) 


0.793(5) 


40 




0.4 


11 




0.321(5) 


0.683(5) 


11 




0.5 


9 




0.617(10) 


0.875(5) 


28 




0.5 


10 




0.439(9) 


0.734(6) 


9.4 




0.5 


11 




0.315(6) 


0.612(6) 


1.6 




0.6 


9 




0.746(11) 


0.893(5) 


20 




0.6 


10 




0.576(9) 


0.767(5) 


19 




0.6 


11 




0.357(6) 


0.584(5) 


55 


m 


0.4 


9 


-0.003(1) 


0.76(4) 


0.40(3) 


8.7 




0.4 


10 


-0.019(4) 


0.47(4) 


0.24(1) 


3.2 




0.4 


11 


0.003(3) 


0.85(6) 


0.45(4) 






0.5 


9 


-0.008(2) 


0.70(4) 


0.42(3) 


29 




0.5 


10 


-0.012(4) 


0.55(5) 


0.31(3) 






0.5 


11 


-0.005(4) 


0.54(5) 


0.28(2) 


0.88 




0.6 


9 


-0.010(2) 


0.69(3) 


0.48(3) 


0.05 




0.6 


10 


-0.006(4) 


0.68(6) 


0.59(6) 






0.6 


11 


-0.055(12) 


0.28(4) 


0.269(4) 





Table 2. Parameters for the j ^ extrapolation of (qq) in the 
transition region, using a linear ((5)), power-law and power 
+ constant Q Ansatz. All four values of j have been used. 
Where there is a star in the column for /Ndf, the central fit 
value was outside the 68% confidence interval, and the quoted 
value is instead taken to be in the middle of the 68% confidence 
interval. 
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Table 3. Inflection points Ts{j) for {qq}{T) at ja = 0.04, 0.02 
and critical temperature Ts obtained from extrapolating Ts{j) 
to j — 0. The uncertainties are estimates of the systematic 
uncertainty in determining the inflection points and in the j — >■ 
extrapolation. 
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Table 4. Estimates for the deconfinement crossover tempera- 
ture Td from the Polyakov loop at ja — 0.04. The fi = result 
is taken from 7 . 
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Fig. 2. The renormahsed Polyakov loop (L) as a function of 
temperature T for fia — 0.35, 0.4, 0.5 and 0.6, in addition to 
the fj, = j = results from T . The open symbols are results 
for ja = 0.04; the shaded symbols are for ja = 0.02. 

3.2 Deconfinement transition 

Figure [2] shows the traditional order parameter for dccon- 
fincmcnt, the Polyakov loop (L), as a function of tempera- 
ture for different values of /i. The Polyakov loop has been 
renormahsed using the same method as described in [71 
I23j . by imposing the condition that L(ar= -j, /i = 0) = 1. 
This determines the renormalisation constant in the 
relation between the bare Polyakov loop Lq and the renor- 
mahsed Polyakov loop L^, 



1 



(5) 



At all /I, we see a transition from a low-temperature 
confined region to a high-temperature deconfined region. 
In contrast to the diquark condensate, we see a clear, sys- 
tematic shift in the transition region towards lower tem- 
peratures as the chemical potential increases. 

For all four /i-values, the Polyakov loop shows a nearly 
linear rise as a function of temperature in a broad region, 
suggesting that the transition is a smooth crossover rather 
than a true phase transition. For this reason, it is diffi- 
cult if not impossible with our present data to identify an 
inflection point which we could use as a definition of a 
pseudocritical deconflnement temperature. We choose in- 
stead to define the crossover region as the region where 



the Polyakov loop at ja — 0.04 rises approximately lin- 
early with T. From Fig. [2] we see that the value of (L) 
increases as j is reduced, and at /la = 0.6 at least, the 
crossover region will most likely move to smaller T in the 
j — > limit. However, we do not have sufhcient statistics 
for ja = 0.02 at low T to make any quantitative statement 
about this. 

In the low-temperature region, an accurate determina- 
tion of the renormahsed Polyakov loop becomes increas- 
ingly difficult because of the exponential growth of the 
renormalisation factor Z^^ . Extending the data in fig. [5] 
to Ta < 0.05 {Nt > 20) is beyond our present capabilities. 

Our summary of transition temperatures taken from 
the ja — 0.04 data is given in table HI 



4 Static quark potential 

The potential between two static quarks (or a quark- 
antiquark pair), and in particular its asymptotic behaviour 
at large separations, has traditionally been taken as the 
tell-tale indicator, or even definition, of confinement of 
quarks [21] • A linearly rising potential has been observed 
in numerous lattice simulations, and has also formed the 
basis of successful phenomenological descriptions of bound 
states of heavy quarks. In QCD with dynamical quarks, 
the string will break at a finite distance, but at interme- 
diate distances a linear rise can still be observed. 

At high temperature, the potential is expected to ex- 
hibit Debye screening, and this has been observed in nu- 
merous calculations of the quark-antiquark free energy us- 
ing Polyakov loop correlators. However, it is not yet clear 
how this quantity relates to the (complex) potential that 
appears in effective theories of heavy quarkonia at high 
temperature [25l[26l[27l[28] . Very recently, the static quark 
potential has also been determined from Wilson loops at 
high temperature [291 : this does not show any screening 
for T < Tc. 

There has also been some recent progress in determin- 
ing the potential between heavy (finite mass) quarks at 
zero [3D] and non-zero |3T] temperature. Some properties 
of bound states of heavy quarks in QC2D at nonzero tem- 
perature and density were reported in |32] : a potential 
model description should reproduce these results. Here we 
compute the static quark potential from Wilson loops for 
our lowest temperature, the 12'^ x 24 lattices. 

In fig. Owe show the static quark potential computed 
from the Wilson loop at Nr = 24, for fia = 0.3, 0.5, 0.7, 0.9 
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Fig. 3. The static quark potential computed from the Wilson 
loop, for the 12^ x 24 lattice and different values of fi, with 
ja = 0.04. 



and 1.1. We find that as we enter the superfluid region, the 
potential becomes slightly flatter, but that this is reversed 
as /X is increased further, leading to a strongly enhanced 
string tension at fia = 0.9, which according to our analysis 
of the Polyakov loop should be in the deconfined region. 
This agrees with the pattern that was already observed in 

0- . . . 

To quantify the variation of the static quark potential 

with /i, we have performed a fit to the Cornell potential. 



V{r)=C{iiJ)+a{fi,j)r + 



(6) 



for each /i and j. The results for the "string tension" a 
and the Coulomb factor a are shown in fig. ID For jia < 0.6 
we do not see any systematic trend, but for larger fi we 
find that a increases rapidly and roughly linearly with /i, 
accompanied by a similar decrease in a. We also see that 
there is no significant dependence on the diquark source 

j- 

The lack of screening at large fi is consistent with what 
was observed at high T in [29 . At present wc do however 
not have a good understanding of why the static quark 
potential should become antiscreened at large /i. There 
are several possibilities: 

1. The transition may be to a medium characterised by 
long-range interactions, rather than by colour screen- 
ing. If that is the case, one would expect typical cor- 
relation lengths to also grow with /u. Reconciling a lin- 
early rising potential with a nonzero Polyakov loop 
remains a challenge, however. 

2. The quark-antiquark potential is screened at large dis- 
tances, but this is not observed in the Wilson loop be- 
cause of poor overlap with the relevant states. This cor- 
responds to the standard scenario at low temperature, 
where explicit mesonic states must be introduced to 
observe string breaking |33j . Clearly, our observation 
of a linearly rising potential (area law for the Wilson 
loop) does not prove that the medium is confining. 
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Fig. 4. The parameters a (top) and a (bottom) in the Cornell 
potential ((B)), as a function of the chemical potential /x, for the 
12^ X 24 lattice. The error bars are statistical uncertainties, 
while the bands indicate systematic uncertainties estimated 
from varying the fitting range in R. 



The increasing slope could however be indicative of a 
large internal energy for static quark-antiquark pairs 
at intermediate distances. 

In ^35^ it was found that the binding energy of heavy 
quarkonia in QC2D increases up to /la ~ 0.7, and de- 
creases again beyond that. This appears to run counter 
to the conjecture above, in which the quarkonium might 
be expected to become more strongly bound. 
3. A more pessimistic scenario is that the whole region 
of > 0.6 could be dominated by lattice artefacts. 
Unfortunately our previous data at (3 = 1.7 are 
probably outside the scaling region so a comparison 
with those is likely to be not very revealing. Simula- 
tions at smaller a, which are underway, should confirm 
or rule out this scenario. 

Computing the static quark potential using Polyakov 
loop correlators rather than Wilson loops might yield fur- 
ther insight into this issue. However, the Polyakov loop 
correlator suffers from the same signal to noise problem 
at low T (large Nr) as the Polyakov loop itself, and we 
have not been able to obtain any signal for Nr = 24 ex- 
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Fig. 5. The gluon dressing function in the vacuum from the 
16^ X 24 and 12'' x 24 lattice. The inset shows the gluon prop- 
agator for infrared momenta. A cylinder cut has been applied 
to the data to reduce lattice spacing artefacts. 



cept for fj,a > 0.9. Results for higher temperatures using 
both the Wilson loop and Polyakov loop correlators will 
be presented in a future publication. 



5 Gluon propagator 

In this section we extend previous studies [51 1^1^ and 
present results for in-medium gluon propagators, where we 
study the dependence on both parameters, chemical po- 
tential and temperature. In Landau gauge only the trans- 
verse part of the vacuum propagator is non-zero. However, 
the external parameters break manifest Lorentz invari- 
ance, hence the gluon propagator D must be decomposed 
into chromoelectric and chromomagnetic modes, De and 
Dm, respectively. 



Dp 



(7) 



where the individual dependence on (discrete) temporal 
and spatial momenta has been omitted. The projectors 
on the longitudinal and transversal spatial subspaces, 
and Pal, are defined by 



Puu{qo,q] 



g2 



(8) 



We have fixed our gauge configurations to Landau gauge 
by maximising the gauge fixing functional 

F[U,g] = J2 Trf/^(x) =^Tr5(x)C/^(x)5t(2;-f A), 

using the standard overrelaxation algorithm. The Lan- 
dau gauge condition has been imposed with a precision 
\d^A^\ < 10-1". 
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Fig. 7. The zero Matsubara mode of the magnetic (top) and 
electric (bottom) gluon propagator as a function of chemical 
potential /i for selected values of the spatial momentum qa , for 
the 16^ X 8 lattice. 



First, in figure [5] we show the gluon propagator and 
dressing function in the vacuum for our two volumes, 12^ x 
24 and 16'^ x 24. Comparing the data for the two volumes, 
we see that finite volume effects are modest for these lat- 
tices. In order to reduce ultraviolet lattice artefacts, we 
have applied a weak cylinder cut [31] . The propagator ex- 
hibits the usual infrared suppression observed in other lat- 
tice studies. 

In figure [5] we show the two lowest Matsubara modes 
for selected spatial momenta as a function of chemical 
potential from the Nr = 24 lattices for different spatial 
volumes. We find at most a very mild volume depen- 
dence, even for zero spatial momentum, with some indi- 
cation that the magnetic propagator is slightly smaller on 
the larger volume. Note that because the available (dis- 
crete) momentum values depend on the spatial volume, 
the selected momenta from the 12'^ and 16'^ lattices do not 
match precisely, and the discrepancy between the propa- 
gator values on the two lattices at nonzero spatial momen- 
tum Qs = \q\ is likely to be at least as much due to the 
slightly different values of Qs as to finite volume effects. 
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Fig. 6. The zeroth (top) and first (bottom) Matsubara mode of the magnetic (left) and electric (right) gluon propagator as a 
function of chemical potential /i for selected values of the spatial momentum qs, for A^^ = 24, different spatial volumes. 



With respect to the infrared suppressed vacuum prop- 
agator shown in fig. [S] we find a mild enhancement at in- 
termediate /i, i.e. in the superfiuid, confined phase, but 
a suppression in the deconfined phase, i.e. for large /i for 
both tensor structures. 

Figure [7] shows the lowest Matsubara mode for the 
high-temperature, 16^ x 8 lattice, again as a function of 
chemical potential and for several different spatial mo- 
menta. Here we find a considerably more complex picture. 
The electric form factor becomes progressively more sup- 
pressed with increasing /x for all momenta, while for the 
magnetic form factor the lowest momentum modes show 
an interesting behaviour, with a peak at /ia « 0.5. For 
large spatial momenta this form factor is instead enhanced 
at large fi. 

We now turn to the thermal behaviour of the gluon 
propagator at fixed chemical potential. Fig. [5] shows the 
zeroth Matsubara modes of the propagators for /xa — 0.5 
and ja = 0.04 on a spatial 16^ lattice as a function of tem- 
perature, with higher modes shown in fig.O The magnetic 
component is slightly enhanced by the effect of chemi- 
cal potential with respect to the vacuum, but interest- 
ingly it hardly feels thermal effects over a large range of 
low to intermediate temperatures, above which it becomes 



screened. The electric propagator is suppressed for lower 
temperatures already. In pure gauge theory [35.36,371 1381 
l39ll40l[4T| the zero mode of the electric gluon shows strong 
enhancement for temperatures below the deconfinement 
transition, whereas the magnetic gluon is screened also be- 
low that transition. However fig. [5] entails that the chosen 
value of n drives the system close to the phase transition 
already at Nt = 20, and hence electric gluon enhance- 
ment may not be observed here. The interplay of tempera- 
ture with chemical potential may also trigger the observed 
non-suppression for the magnetic mode. Our results are 
however in qualitative agreement with those obtained in 
a recent study of QCD with twisted-mass fermions |42) . 
A more detailed analysis of both medium effects and in 
particular their mutual interaction will be presented else- 
where. 



To further quantify the variation of the gluon prop- 
agator with T and /x we employ a global fit D^j^^lq^) 

to all modes, with = + Qq- The functional form we 
use is inspired by [35] ■ However, for the momenta at hand 
the perturbative running of [36J can be neglected. In the 



8 



Tamer Boz et al.: Phase transitions in 2-colour matter 



■ 




= 24,n = 0=j 


< 




= 24 


A 




= 16 


O 




= 12 


□ 




= 8 


o 




= 4 








= 24,n = 0=j 


< 




= 24 


A 




= 16 


O 




= 12 


□ 




= 8 


o 


N 

T 


= 4 




Thermal behaviour of the zeroth Matsubara mode of 
netic (top) and electric (bottom) propagators at — 



Fig. 8. 

the ma] 

0.5 and ja = 0.04 on 16^ x Nr lattices 



54- 















s 


A 






W 


A * 
X 




A 

4 


W 



I 

0.5 




+ 


0* 


mode, 


= 24 


V 


l" 


mode, N 


= 24 


A 


l" 


mode, N 


= 16 


< 


l" 


mode, N 

t 


= 8 




2«' 


mode, N 


= 24 


A 




mode, N 


= 16 




2"' 


mode. 


= 8 



-5 

1 


I 






35 
D 


X 


■ 




a X 
O a 


CD 


■ 
@ 


o 









X 


0'^ 


mode, N = 

X 


24 




O 


l" 


mode, N = 


24 




□ 


l" 


mode, N = 


16 




o 


l" 


mode, = 


8 




o 


2°" 


mode, N = 


24 




n 


2°^ 


mode, N = 


16 




® 


2" 


mode, N = 


8 





0|g 

eg 



1.5 



Fig. 9. Thermal behaviour of higher Matsubara modes of the 
magnetic (top) and electric (bottom) gluon propagators. 



vacuum this gives us a three-parameter fit 

D^Eil') - J^T^il'+A'-M/E)-'"'- . (10) 

At = = j we find Aa = 0.999(3), um = as = 6.85(3) 
and &M = &i; = -1-031(2) on the 16^ x 24 lattice, with a 
per degree of freedom of around 8 for the magnetic mode 
and 5 for the exectric mode. There is a slight volume de- 
pendence, with the 12^ x 24 lattice yielding A = 0.961(5). 
The normalisation A is taken to be independent of T and 
(U, but medium effects modify Um/e and 6m/_e for mag- 
netic and electric modes individually. The results for the 
fit parameters are shown in fig. [TU] as functions of A',- on 
the 16^ X Nr lattices at ^a = 0.5 and ja = 0.04, and in fig. 
[TT] as functions of and lattice volume. For the available 
data we have found the dependence on j to be weak. At 
A^^ = 5 we did not obtain any satisfactory fit for the mag- 
netic form factor, so these points are absent from fig. [TUl 
To illustrate the quality of these fits, fig. fT2l shows the 
fits for fia = 0.5, ja = 0.04 on the 16^ x 24 lattice. We see 
that ([TO)) gives a reasonable description for both modes, 



but the underlying Ansatz that the gluon propagator is 
a function of the four-momentum only, works less well 
for the magnetic form factor. The /^df is around 10 for 
most fits, with the electric form factor giving in general a 
somewhat better x^- We are also planning to employ fit 
models inspired by hard dense loop perturbation theory, 
which depends separately on and qg. 



6 Conclusions 

We have carried out a detailed investigation of several 
aspects of 2-colour, 2-flavour QCD with m-^/mp « 0.8 at 
nonzero temperature T and quark chemical potential ^. 
Our main findings are summarised below. 

1. We have located the superfluid to normal and decon- 
finement transitions in the region 0.35 < /la < 0.6 
(/J, = 385 — 665 MeV). The superfluid to normal tran- 
sition temperature is remarkably constant in this 
region, while the deconfinement temperature Td shows 
a decrease with /i which appears to continuously con- 
nect to the ^ = Q transition identified in [7]. It also 




Fig. 10. A'^^ dependence of fit parameters aj^i/^ and fcjvi/B for 
16^ X Nt lattices at /xa = 0.5 and ja — 0.04. 



appears to extrapolate smoothly to the high-/i, low-T 
transition previously observed OIHIII] , although in the 
absence on any accurate data for the Polyakov loop at 
low temperature this must be taken merely as indica- 
tive. 

2. The superfluid to normal transition appears to behave 
like a second order phase transition, while the decon- 
finement transition looks like a smooth crossover. This 
would have to be backed up with a careful finite vol- 
ume and critical scaling analysis. 

3. The static quark potential at low temperature appears 
to become more, rather than less confining at large /i, 
suggesting that the dense medium with (L) 7^ is not 
an ordinary, deconfined quark-gluon plasma. It must 
however be noted that string breaking will occur even 
though it is not seen in the potential obtained from 
the Wilson loop because of the poor overlap with the 
relevant states. 

4. The electric (longitudinal) gluon propagator in Landau 
gauge becomes strongly screened with increasing tem- 
perature and chemical potential. The magnetic (trans- 
verse) gluon shows little sensitivity to temperature, 
and exhibits a mild enhancement at intermediate /i 
before coming suppressed at large 

The structure of the phase diagram is summarised in 
fig- HSl We also include the estimates for the deconfine- 
ment transition given in [7]. We see that these are con- 
sistent with, but slightly above, our new estimates in this 
paper. The indications are that the deconfinement line 
crosses the superfluid to normal transition line, giving rise 
to a region of deconfined, superfluid matter, but in the ab- 
sence of precise data at larger /i, we can not say this with 
any degree of certainty. To map out the phase boundaries 
with greater precision and to clearly establish the order of 
the transitions, large-scale simulations on several spatial 
volumes will be required. This goes beyond our current 
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Fig. 11. Dependence of fit parameters aM/E (top) and bM/E 
(bottom) on and lattice volume. 



computational capabilities, but is an interesting topic for 
future investigations. 

The lack of screening (and even antiscreening) observed 
in the static quark potential requires further investigation 
to establish whether this is a signal of an exotic state of 
matter or a result of poor overlap with relevant states 
and/or lattice artefacts. We are planning to compute the 
static quark potential at higher temperatures, using both 
the Wilson loop and Polyakov loop correlator, to further 
elucidate this. 

The screening of the static magnetic gluon propagator 
(zero Matsubara mode) at high /i and low T is a clear 
signal of the breakdown of resummed perturbation theory, 
which predicts that this mode is unscreened to all orders. 
This is likely to be relevant also for real QCD, and should 
be taken into account in any future study using model 
gluon propagators to compute for example the superfluid 
or superconducting gap. 

In an ongoing study we will compare the lattice gluon 
propagator at non- vanishing chemical potential and tem- 
perature with results from functional continuum methods, 
extending studies of thermal propagators p51HilH5] to fi- 
nite density. As there are no insurmountable problems in 
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aj = 0.04 on the 16^ x 24 lattice. Note that the functions are 
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QC2D, in particular no sign problem, a direct analysis 
serves to identify possible technical limitations in either 
method, stemming from finite size and finite volume arte- 
facts in the lattice formulation, or from inevitably neces- 
sary truncations in the continuum description. We are also 
in the process of computing the quark propagator, which 
will give further input to these studies. 

We are currently extending our study of QC2D to 
smaller lattice spacings, which will enable us to perform a 
controlled extrapolation to the continuum limit and clar- 
ify the possible role of lattice artefacts at large /i. 
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